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Abstract 

A detailed quantitative analysis of the transition process mediated by a sphaleron type 
non-Abelian gauge field configuration in a static Einstein universe is carried out. By 
examining spectra of the fluctuation operators and applying the zeta function regular- 
ization scheme, a closed analytical expression for the transition rate at the one-loop level 
is derived. This is a unique example of an exact solution for a sphaleron model in 3 + 1 
spacetime dimensions. 



1 Introduction 



The discovery of vacuum periodicity |], [0 and sphalerons || in the standard model 
has revealed a considerable fermion number violation in the theory at high tempera- 
tures. The fermion number changes in the processes of barrier transitions between the 
distinct topological sectors of the theory, where the barrier height is determined by the 
sphaleron energy Such processes become unsuppressed at temperatures of the order of 
the sphaleron mass. Since the different topological sectors are in thermal equilibrium 
at high temperatures, the sphaleron transitions lead to the dissipation of the baryon 
asymmetry produced at earlier times, that is, during the electroweak phase transition or 
at GUT energies ||]. The estimates for the transition rate @-||] show that the sphaleron 
mechanism can be efficient enough to reduce the asymmetry by several orders of magni- 
tude. 

The sphaleron solution in the standard model is known only numerically, which 
causes considerable computational difficulties J7J-H]. Other field-theoretical sphaleron 
models have therefore become important |I0|| , [ pH , |12| . Some of these models are exactly 



solvable in the sense that the corresponding sphaleron transition rate can be evaluated 
analytically at the one- loop level IpTf , ||12|| , which provides a closer insight into the physics. 
Unfortunately, all of these models exist only in low spacetime dimensions. There is also 
an alternative example of four dimensional sphalerons which appear in the theory of 
a self-gravitating non-Abelian gauge field. Gravity violates the scale invariance and 
plays therefore a role similar to that of a Higgs field. This results in the existence of 
classical |13j sphaleron-like solutions in the theory [f[4 , fT5| , which are also known only 
numerically. 

The purpose of the present paper is to investigate a sphaleron model that is exactly 
solvable and yet exists in 3 + 1 spacetime dimensions, which distinguishes it from the 
other known models. Similarly to the example mentioned above, this model deals with a 
non-Abelian gauge field interacting with gravity, but now gravity is regarded as a fixed 
external field. Specifically, we shall consider the theory of a pure SU (2) Yang-Mills 
field in a static Einstein universe. This theory admits an analytically known sphaleron 
solution which we shall call cosmological sphaleron. The analysis of the corresponding 
transition problem is the subject of this paper. 

In fact, this type of sphaleron was found long ago. The solution itself has been 
discussed in various contexts, together with the other related solutions of the (Einstein)- 
Yang-Mills field equations |j[6|l-|[LSj|. The sphaleron nature of the solution was also recog- 
nized and discussed in connection with the finite volume QCD and in a cosmological 



context as well |P0|| , However, an analysis of the corresponding sphaleron transition 
problem has been lacking so far. 

The cosmological sphaleron is distinguished by the important property that the 
sphaleron configuration consists of the gauge field alone. This implies that the sphaleron 
is not plagued with the various symmetry breaking phase transitions, and this also 
ensures that the dynamics of the sphaleron mediated processes is conformally invariant 
- up to the anomalous scale dependence of the gauge coupling constant. Since the 
sphaleron has very high symmetry, there is an opportunity to analyze the sphaleron 
transition problem analytically at the one-loop level, which is of great methodological 
interest. (We shall use the names "cosmological" and "universe", although the model 
can be applied not only in the cosmological context.) 
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The problem which will be investigated below can be formulated as follows. Con- 
sider a gauge field in a static Einstein universe at finite temperature. Specifically, consider 
the thermal ensemble over one of the topological vacua of the field. Find the rate of the 
decay of this thermal state due to the diffusion into the neighboring topological sector. 
In order to obtain the answer, we use the Langer-Affleck formula and take only the 
bosonic degrees of freedom into account (the fermion contribution can be considered in 
a similar way). We do not assume the high temperature limit and take the sum over all 
Matsubara modes. We use the zeta function regularization scheme, which allows us to 
entirely carry out the analysis. Our principal results are given by Eqs. (|6.20 )-(S.27) and 
Eq.( |6.37|) and presented in Figs. 1-3. The rest of the paper is organized as follows. The 
basic properties of the model, such as topological vacua and the sphaleron solution are 
discussed in Sec. II. A brief derivation, based on the path integral methods, of the decay 
rate of an unstable phase at finite temperature is given in Sec. III. The path integration 
procedure is outlined in Sec. IV. The spectra of the fluctuation operators are analyzed in 
Sec.V. The determinants of these operators are calculated in Sec. VI which also includes 
the consideration of the high temperature limit. Sec. VII. contains some concluding re- 
marks, and the derivation of numerous formulae used in the zeta function approach is 
given in the Appendix. 

Throughout the paper the units h = c = k b = 1 are used. The symbol g stands 
for the gauge coupling constant, whereas the spacetime metric is denoted by g. 



2 The sphaleron on S 3 

Consider the static Einstein universe (M, g), where M = R 1 x S 3 , and the metric is 

ds 2 = a?(-d V 2 + dnl). (2.1) 

Here a is a constant scale factor, and the line element on S 3 is parameterized by 

<m\ = d? + sin 2 i(d§ 2 + sin 2 §d^ 2 ) , (2.2) 

where £ G [0, n], and (p are the usual spherical coordinates on S 2 . 
The model under consideration is defined by the action 

S[A] = ~tT [ V=gd 4 x, (2.3) 

where A = A^dx^ 1 = TpA^dx 11 is the gauge field, and F^ v = d^A u — d u A^ — i[A^, A u ] is 
the field tensor. The hermitian group generators are T p = t p /2 with r a being the Pauli 
matrices. The metric g in Eq.( [2.3|) is given by Eq.( |2.1|) , such that the action defines 
the theory of a pure non-Abelian SU(2) gauge field in the static Einstein universe. 
(Throughout the paper we neglect the back reaction of the gauge field on the spacetime 
geometry, which can be justified if a ^> l p i/g, where l p i is Planck's length.) The classical 
equations of motion following from the action are 

-t[A fl ,F^}=0, (2.4) 

where V M is the covariant derivative with respect to the spacetime metric. 
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First, we need to describe the topological vacua of the gauge field in this case. We do 
so by introducing a smooth, time-independent function on the manifold, U(x) e SU(2), 
where x = x m , m = 1, 2, 3, thus defining the mapping 

*7(x) : S 3 -> SU{2). (2.5) 

Any such mapping can be characterized by an integer winding number 

k[U] = — L tr [ UdU- 1 A UdU~ l A UdU' 1 , (2.6) 

247T JS 3 

such that the set of all U's falls into a countable sequence of disjoint homotopy classes. 
The representative of the k-th class, U^ k \ k[Z7^] = k, can be chosen as 

I/W(x) = U {k \£,$,<p) = exp{-ik£ n a r a }, (2.7) 

where n a = (sin $ cos if, sin t? sin ip, cos i?). The functions C/ generate static gauge trans- 
formations, 

A —> UAU~ l + iUdll^ 1 . (2.8) 

The elements of the zeroth homotopy class, U^°\ give rise to small gauge transformations 
which can be continuously deformed to identity; other functions, U^ k \ k 0, generate 
large transformations. A vacuum of the gauge field is a pure gauge, A vac = iUdU^ 1 . 
Since all U's split into homotopy classes, all pure gauges decompose into disjoint sets 
called topological vacua. The k-th topological vacuum in the temporal gauge is 

A (fe) (x) = iU {k) dU {k) - 1 . (2.9) 

By construction, the Chern-Simons number of this field configuration coincides with the 
winding number k. 

Distinct topological vacua can not be joined by a continuous interpolating sequence 
of pure gauge configurations, iUdU' 1 , since this would require a change in the winding 
number k[Z7]. However, one can join them by a family of non- vacuum fields. That is how 
one can see that the model admits a sphaleron solution. Consider the two neighboring 
vacua given by Eqs. (|2.7|) , ( |2.9D : A^ = and A^\ They can by joined by the following 
path in the configuration space: 

A[h]=l l±^uW du m-\ (2.io) 

where the parameter h G [—1,1]. (Applying a large gauge transformation one can reduce 
such a path to a non contractible loop.) The energy-momentum tensor for this field is 

( h 2 - l ) 2 

such that the energy is given by 



T?(A[h}) = \ 2 I diag(3, -1, -1, -1), (2.11) 



E[h] = I T^d 3 x = ^(h 2 - l) 2 . (2.12) 
J g a 

This function has the typical barrier shape: it vanishes at the vacuum values of h, 
h — ±1, and reaches its maximum in between, at h — 0. The top of the barrier relates 
to the field configuration 

=A[h = 0}= i - U {l) dU {1) -\ (2.13) 
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with the energy E max = 3n 2 /g 2 a.. 

Similarly, one can define E max for any other interpolating path, and then minimize 
the result over all paths. If a non zero minimum exists then it relates to an unstable 
classical solution called sphaleron. By construction, the sphaleron energy defines the 



minimal height of the potential barrier H], [10 



To carry out such a program would, however, be too difficult a task which has never 
been done in reality. Instead, we simply check that the field fl2.13|) solves the classical 



equations of motion. To see this, it is illuminating to allow the parameter h in Eq. ( |2.10|) 



to depend on time, h — > h(rj). Then the Yang-Mills equations ( j2.4|) for the field Q2.1U ) 
reduce to one non-trivial equation, which admits a first integral: 

• 2h(h 2 1) 0. =► (^V + (tf _ I)' = e , (2.14) 



drf 



with e being an integration constant. Effectively, these equations describe a particle 
moving in the one- dimensional double- well potential. When e = 1, one finds the static 
solution h(rj) = 0, which describes an unstable equilibrium of the particle on the top 
of the barrier. This shows that field configuration ( |2.13| ) indeed solves the equations 
of motion and relates to a saddle point of the energy functional, such that it can be 
naturally called sphaleron. Later we shall see that this solution has only one unstable 
mode. In addition, for e = 1, Eqs. Q2.14l) admit the solution for the particle rolling 
down the barrier, h{rj) = \/2/ cosh v2(t? — Vo), which describes the time evolution of the 
sphaleron during its classical decay. 

Of course, these arguments do not prove that the sphaleron relates to the absolute 
minimum of energy for static, nonvacuum solutions. Notice however, that ( |2.13|) is the 



only static, nonvacuum, 5*0(4) symmetric solution | 17fl . It is therefore very plausible 
that this solution does indeed minimize the energy. 

Another handy form for the sphaleron solution (|2.13 ) can be achieved as follows. 
Introduce the left and right invariant 1-forms on S 3 , 

uj a L = '-tr^U^dU^- 1 ), u a R = ^tr(r a U m - 1 dU m ), (2.15) 
which satisfy the Maurer-Cartan equations 

du a + e abc u h A uj c = 0. (2.16) 
This allows us to represent the field ( |2.13| ) as 

A (sp) =T a u a L . (2.17) 



It is worth noting the following feature of this solution: the sphaleron configura- 
tion consists of the gauge field alone. This, together with the high symmetry of the 
solution, will be of crucial importance for the analysis below. It is well known that in 
the Minkowski space, the existence of static, finite energy solutions for the pure gauge 
field is ruled out by the scaling arguments. However, these arguments do not generally 
apply in curved spacetime, since the invariance with respect to the rescaling of the co- 
ordinates, x — > Ax, is broken by the curvature. It is therefore gravity which ensures the 
existence of the static sphaleron solution. Since the spacetime geometry is homogeneous 



4 



and isotropic, the sphaleron inherits the same symmetries, such that, for instance, its 
energy-momentum tensor has the manifest S'0(4)-symmetric structure. 

Let us mention also the following point: since the Yang- Mills equations are confor- 
mally invariant, and the geometry (|2.1| ) is conformally flat, there exists the Minkowski 
space counterpart of the sphaleron solution. However, this flat spacetime solution is, 
of course, no longer static. Specifically, if one chooses the conformal factor of the 
metric ( |2.1| ) as a = a(r/, £) = (cos 77 + cos£) _1 and introduces the new coordinates 
£±r = tan((?7±£)/2), the metric assumes the standard flat form [E0|| . In the new coordi- 



nates, the sphaleron field (|2.13| ) becomes a member of the family of the elliptic solutions 



22], which describe spherical shells of the Yang-Mills radiation in the Minkowski space 
23 1 - In this sense, one can think of the sphaleron as being a radiative solution which is 
rendered static by gravity. 

We will also need some knowledge about instantons in the model. We first note 
that the vacua and the sphaleron, since they are static, can be regarded as Euclidean 
solutions. Next, let us pass to the imaginary time in Eq. fl2.14p , rj — > —ir, 

= -e. (2.18) 

Apart from the static solutions, this equation admits also the interpolating solutions for 
e = 0, such that h(—oo) = —1, h(oo) = 1. The periodic solutions exist for < e < 1, 
and the corresponding period is bounded from below, 

t > V2n. (2.19) 

Other known instanton solutions on S 3 can be obtained from the flat space BPST in- 




stantons by making use of the conformal invariance of the YM equations [19 . 



3 The sphaleron transition rate 

Consider the low energy excitations over the k-th topological vacuum, A^ k '('x) — > 
A«(x,t) = A«(x) + <L4(x,t). If the energy is small compared to the barrier height, 
3ir 2 /g 2 a, then the excitations over the distinct vacua are classically independent. There 
is a nonzero amplitude for the quantum tunneling between distinct sectors, however the 
corresponding probability is exponentially small. On the perturbative level, one can 
consider the excitations in each sector independently. The energy of the ground state 
excitation in each sector is 1/a (see Eq. (|5.20Q below). The necessary condition for the 
smallness of the energy of the excitations is therefore g 2 /37r 2 <C 1. 

Consider the zeroth topological sector and assume a thermal distribution for the 
states in this sector. There is a finite probability for such a thermal system to decay, 
both because of the underbarrier tunneling and due to the overbarrier thermal excitation. 
According to the Langer- Affleck theory of the metastable phase [ESI, the decay rate is 



proportional to the imaginary part of the free energy. To estimate the latter, it is 
convenient to use the path integral approach (the precise definition of the path integration 
procedure will be given in the next section). 

The partition function of the gauge field is 



Z = exp(-/3F) = J d[A] exp(-S E [A]). (3.1) 
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In this expression, <Se[A] is the Euclidean action of the gauge field in the static Rieman- 
nian space (S 1 x S 3 , g), with g being the analytic continuation of the metric (2.1) to the 



imaginary time, r] — > —it, t G [0, 0\. In the weak coupling limit, one can approximate 
the partition function by the sum over the classical extrema, as Z ~ J2j Zj, where 
the semiclassical contribution of the j-th extremum is 

Zj = exp{-PFj) = exp(-S[A^ j} }) J d[<p] exp(-5 2 S' j ). (3.2) 

The action for the fluctuations around the j-th extremum, A^ — > A^ + (p, can be 
represented as 

S[A& + <p]= S[A&] + 5 2 S 3 + . . . , 5 2 S 3 = f dr\ (<p, D^) ^/gd 3 x, (3.3) 

Jo J s 3 

where D-,- = DfA^] is the Gaussian fluctuation operator. This gives the one loop 
expression for the partition function, 

ZkEZj = E !5£24»!). (3 . 4) 

Assume that this sum is dominated by two terms, Z ~ Zq + Zi, where Zq and Z\ are the 
contributions of the vacuum and the sphaleron, respectively. Other periodic instantons 
that could exist for a given value of j3 are assumed to have a large action. 

The sphaleron fluctuation operator Di has at least one negative eigenvalue uPi < 0. 
Under the condition specified by the lower bound in Eq. (|3.8| ) below, there is only one 
negative eigenvalue. This implies that Z\ is purely imaginary, and the free energy of the 
whole system picks up the imaginary part 

Im(/3F) = -Im hiZ ~ -Imln (z (l + —)) ~ -— Im Z x . (3.5) 

V Zq J Zq 



According to Langer |Pq|,the imaginary part of the free energy is to be interpreted 
as giving rise to the decay rate of the unstable phase built over the perturbative vacuum 
as 

r = M Im F, (3.6) 

where the damping constant k is the real time decay rate of the sphaleron configuration in 
the heat bath. In the weak coupling limit one has \k\ = |u;_| 0, which finally determines 
the decay rate to be 

_. |a?_| ImZi 



7T Zq 

This formula holds in the following range of temperatures EE 



(3.7) 



\u-\ 1 3tt 2 

< S « -JT- 3. 



2vr (3 g 

The lower bound rules out the periodic instantons which play the leading role at low 
temperatures p8 (for the cosmo logical sphaleron one has \lu_\ = y/2, such that this con- 



dition is opposite to that specified by fl2.19| )). The upper bound is the sphaleron energy 
which must exceed the temperature, otherwise the system would not be metastable. No- 
tice that this formula involves only conformally invariant quantities, the conformal factor 
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a drops out. Another crucial assumption is the weak coupling limit, g 2 /4:7i <C 1. First of 
all, this ensures the very existence of the thermal ensemble. In addition, it justifies the 
validity of the Gaussian approximation, and, moreover, leads to the weak damping limit 
for k. 



4 The path integration procedure 

We outline below the main steps of the path integration procedure. It is worth noting 
that the gauge field theory on S 3 resembles that on S A . We shall therefore mainly follow 
the approach given in Refs. ||25|| — 127|. 



Passing to the imaginary time, the spacetime metric ( |2.1| ) becomes 

ds 2 = a 2 (rfr 2 + dQ 2 3 ), (4.1) 

where r G [0,(3]. We assume that the coordinates are dimensionless, implying that 
[a] = [L\. Notice that 1/(3 is the conformal temperature. The physical temperature T is 
defined with respect to the physical time, ar, such that T = l/(3a. 
The Euclidean action of the gauge field is 

S E [A] = ±tr J ^gd 4 x > 0. (4.2) 

Consider small fluctuations around the j-th. extremum of the action — > A^p + cp^, 
where the values j = 1,0 refer to the sphaleron and the vacuum configurations, respec- 
tively (we shall omit this index where possible). The infinitesimal gauge transformations 
act as y? M — > ip^f* = (p^ + D^a, where D^a = V^a — i[A^ a], and a is a Lie algebra 
valued scalar field. Define the following operators: 

Y)ip u = JsA(p u + D v [D a ^f\ M<p v = -D a D a ip v + Rlip u + 2i[F v „<p a \. (4.3) 

These are the vector fluctuation operator and the gauge fixed fluctuation operator, re- 
spectively. Introduce also the Faddeev-Popov operator 

M FP a = -D a D a a. (4.4) 

In these formulas R u a is the Ricci tensor for the geometry ( |4.1| ) , and F va is the background 
gauge field tensor. These operators are self-adjoint (symmetric) with respect to the scalar 
products 

(<p, <p') = 2tr J (p a cp' a ^/gd 4 x, (a, a') = 2 tr J aa'y/gd 4 x. (4.5) 



The norms are \\ip\\ = yj (ip, ip), and ||a|| = yj (a, a). The action can be expanded 



as 



S s [A + <p]K S E [A] + ±(<p,f><p). (4.6) 



The Gaussian path integral is 

Z = exp(-S E [A}) J ^>]exp (-^(^f)^ , (4.7) 
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where the Faddeev-Popov measure is 

d FP [<p] = d[<p]g{<p)F(<p), 1 = Q(<p) J d[a}Hv {a) ), (4.8) 
and T is the gauge fixing function. The fluctuations if can be decomposed as 

tp„ = D^a + (4.9) 

where the pure gauge part D^a is annihilated by D, and £ M is orthogonal to all gauge 
modes, D a £ > (T = 0. The fields £ M and a can be expanded with respect to the eigenf unctions 
of D and M FP : 



5> 



a. 



r>C k = A fc ejfe, M FP a n = q n a 



The gauge fixing function is chosen to be 

1 



J 7 (ip) = exp 



2g'- 



1 



2r 



(4.10) 



(4.11) 



and the integration measure is the square root of the determinant of the metric on the 

- dB n y/ q n I \ a n \ I , 



function space [27]: 

%]=n^^ii^iiri / ' n 



d[a] = l[-^dB n \\a n \\. (4.12) 



2ng 



2ng 



2ng 



Here jjq/ {y2wg) is a normalization factor with being an arbitrary normalization scale; 
the prime indicates that terms with q n = should be omitted. Taking (|4.8|) - (|4.12|) into 
account, the Gaussian path integral in (PO) reduces to 



(4.13) 



J k V2irg n V2ng { 2g \ k 



where ^Jq^ has been absorbed in B n , and 



<r = fn 



/'o 



27rg 



dB n \\a n \ | exp 



2g'- 



J2 B nQn\\ a n\ 



(4.14) 



Let us first apply these formulae to the sphaleron. In this case, as we will see in the 
next section, the vector fluctuation operator has one negative eigenvalue, A_ < 0, whereas 
all the other eigenvalues Xf., q n in (|4-13|) , (|4.14|) are positive (zero modes are absent). The 
integral over in ( (4.13| ) can be defined by analytic continuation [2"3|, and the result is 
/io/2z^/|A_|. The rest are well-defined Gaussian integrals. It is convenient to introduce 

the conformally invariant dimensionless operators M and M FP whose eigenvalues are io\ 
and Lot: 



M = i-M, M FP = \m fp ; 



lot 



At. — — -, 



q =^ 

Hn 9 • 



(4.15) 



This implies that the partition function for the Gaussian fluctuations around the 
sphaleron is 



2i\j 
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where Det' has all nonpositive eigenvalues omitted, and the index 1 referring to the 
sphaleron is restored. Notice that Det'(M) must be computed on the space of all vector 
fluctuations ip", and not only for those satisfying D a ip a = 0. 

Consider now the vacuum case. As we will see below, the ghost operator M FP 
in this case has three zero modes a p = t p /2 (p = 1,2,3) related to the global gauge 
rotations of the vacuum A = 0. The norm of these modes is 

\\a p \\ 2 =[ drf y/gd 3 x = 2vr 2 /?a 4 . (4.17) 
Jo J s s 

The quantity Q specified by Q4.14 ) then reads 

^ n ft f \ ft 



with 

T = f[ / ^dB p \\a p \\ = V SU(2) ^^P^, (4.19) 



where the integration over B p gives the volume of the stability group, in our normalization 
it is Vsu(2) — 167T 2 . In addition, there are three constant vector modes annihilated 
by the vector fluctuation operator: = 5qT p /2, \\(, p \\ 2 = 27r 2 /3a 2 . Under the gauge 
transformation generated by 

U p (r) = exp (ir^-hn , U p (0) = U P (J3), I G Z, (4.20) 

these modes change according to 

U P (r) : e p = - (l + — J . (4.21) 

This shows that the range of integration over dC p in ( [4.13 ) is finite, C p G [0, 4ff//3], The 
contribution of these three modes to (|4.13 ) therefore is 



All other eigenvalues are positive, which finally gives for the fluctuations around the 
vacuum 

4tt Det'{M FP //i 2 a 2 ) 

We shall omit below the factor /x , ft a ~^ a ; such that a will be understood as the radius 
of the universe expressed in units of an arbitrary length scale. 



Z ° ~ ,3 a 3fl3 / - - = - • ^- I6 > 



5 Spectra of the fluctuation operators 

To analyze the spectra of the conformally invariant operators M and M FP defined by 
Eqs. (|4.3| ), (|4.4|) , (|4.15 ), one can put a = 1 in the line element. We introduce the 1-form 
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basis u a } on the spacetime manifold, where u° = dr, and uj a = io a L are the left 
invariant 1-forms given by Eq.( [2.15 ). The metric is 

d s 2 = u° ® u° + cu a <g> uj a . (5.1) 

Let {eo, e a } be the corresponding dual tetrad; here e a = e% are the left invariant vector 
fields on S 3 . Introduce the right invariant fields dual to the 1-forms cjfj. The following 
commutation relations hold 

[ e ai e b] = 2 £ abc£ci [ e ai e f ] = ^ £ abc^c ', [ e a ) ^ = 0- (5-2) 

Let Vo and V a denote the covariant derivatives along the tetrad vectors {eo,e a }. The 
following tetrad rotation coefficients do not vanish: 

V a eft = Sabc^ci V a u/ = S abc LO c . (5.3) 

Let us represent the gauge field of the vacuum and the sphaleron as 

A^=ju s T s , (5.4) 



where j = 0, 1, respectively. All this suggests to expand the fluctuations as JTT) 

<p = (J)°uj + <f> a u a , (5.5) 

where 0° and <p a are the scalar and the vector fluctuations, respectively. 

Using Eqs.( [4.3| ), (|4.4|) , ( |5.1| )- (|5.5| ), a straightforward calculation shows that the 
fluctuation operators Mj (j = 0, 1) decompose into the direct sum of the two operators 
acting on the scalar and the vector fluctuations, respectively: 

M j( j) = (Mf a V°) uj° + (M] ec <p a ) oo a . (5.6) 

Here the scalar fluctuation operator M scal formally coincides with the ghost operator 
M FP introduced in the previous section, 

j^coy = _( Vo Vo + V a V a )0° - j {2V S [T S , <f>°] + [T s , [T s , 0°]]} . (5.7) 

The vector operator reads 

Mf c $ a = -(V Vo + V C V C - 4)0 a - 2e abc V b <p c + 

+ j {-2V S [T S , <j) a } + Ae abs [T s , <j> b } - [T„ [T s , a ]]} . (5.8) 
Each of these operators decomposes into the direct sum of a temporal and a spatial part, 

M = -^- 2 +M, (5.9) 



such that the problem reduces to the study of the corresponding spatial operators J\4^ cal 
Cl vec 

3 ' 

We now introduce 



and Mf c 



L a = Z -e^, t a = l -e^ L 2 = L a L a = L a L a , (5.10) 
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which are the SO (4) angular momentum operators, since 

[L a , = ie a bcL c , [L a , L b ] = 2£ afec L c , [L a ,L&]=0. (5-H) 

The commuting operators are L 2 , L3 and L3. The eigenvalues are similar to those for the 
SO (3) case, but the angular momentum can now assume both integer and half-integer 
values: 

13 

L 2 = /(/ + l); I — 0,-,l, L 3 = m, L 3 = rh; m, m — — I, — I + 1, . . . I. (5.12) 

Next, expanding the fluctuations over the basis of the Lie algebra, 4>° = 4>pT P , <p a = 4>pT P , 
we define the spin and isospin operators S and T by 

S a 0j = ~£ a6c </£, T p 0? = T p 0° = ~ £pr X (5.13) 

which satisfy the usual commutation relations. One has S 2 = T 2 = 2, which corresponds 
to the unit spin and unit isospin, respectively. 

Using Eqs.(|5.7|)-( |5.13j ) one can represent the spatial operators for the fluctuations 
around the sphaleron (j = 1) and the vacuum (j = 0) configurations as 

MT = 2 (f? + (L + S + f ) 2 - l) , M{ cal = 2 (f? + (L + f ) 2 - l) , 

Ml ec = 2 (L 2 + (L + S) 2 ) , Ml cal = At 2 . (5.14) 

The fluctuation operators therefore reduce to the combinations of the angular momentum 
operators whose spectra can be analyzed by the usual methods. Let us illustrate the 
procedure for the operator J\4\ ec . 

We represent the operator in the form Xi\ ec = 2 (h 2 + J 2 — l), where J = L + K, 

and K = L + T. The commuting operators are L 2 , J 2 , L 3 and J3, such that the 
eigenvalues of M.\ ec and their degeneracies read u 2 = 2{/(Z + 1) + j(j + 1) — 1} and 
d = u(2l + 1) (2j + 1), respectively. Here v is the degeneracy factor associated with the 
several possibilities to obtain a given value of j for a given value of /. 

To find let us fix a value I >2. K is the sum of the two unit angular momenta, 
such that K = 0,1, 2. For K = 2 the possible values of j are I — 2,1 — 1 ...,/ + 2, for 
K = 1 one obtains j = I — 1, 1, 1 + 1 and for K = one has j = I. The values j — I ± 2 
therefore appear only when K = 2, the values j — I ± 1 are possible when K = 1 or 
K = 2, and there are three different ways to get j = I, that is, when K = 0,1, 2. Thus, 
if we write j — I + a, where a = 0, ±1, ±2, the degeneracies v are as follows: v = 3 for 
a = 0, v = 2 for a = ±1 and v = 1 for a = ±2. As a result, the eigenvalues of Ai\ ec and 
their degeneracies can be represented as 

uj 2 = (2/+(T+l) 2 +(x 2 -3, d= (3-H)(2M-l)(2Z+2<7+l), I > 2, a = 0,±l,±2. (5.15) 

For / < 2, the factor v changes. For instance for / = one has j = K, such 
that v = 1 for a = 2, 1, and v = for a = —1,-2. Similarly for / = 1/2 one 
obtains v = {1,2,2,0,0} for o = {2,1,0,-1,-2}, respectively. For / = 1 one has 
v = {1, 2, 3, 1, 0} and for I = 3/2 the result is v = {1, 2, 3, 2, 0}. 
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Introducing the new quantum number, n = 21 + a + 1, the whole spectrum can be 
given in the following compact form: 



MT 



2 i 2 

n + a 



d 



a 



a = 0,±1,±2; 



(5.16) 



where v — 3— |er| for n > 3, and v 
operators one similarly obtains 



^<W + 5n2<5iH f° r n — 1,2. For the other fluctuation 



Ml 
M v e 



seal 



10 



2 i 2 

n + a 



2 i 2 



d 



n 



a 



a 



0,±1, n > 2: 



M 



seal 



0J~ 



d = 3(n 2 -a 2 ), cx = 0,±l, n > 2; 
n 2 — 1, d = 3n 2 , n > 1. 



(5.17) 



All these eigenvalues are positive for n > 2. When n — 1, the operator A^i ec has 
one negative mode with eigenvalue uj 2 _ = —2. It is worth noting that the sphaleron 
does not possess any zero modes. This can be understood as follows: the sphaleron 
solution completely shares the symmetries of the 3-space - both are SO (A) symmetric. 
The sphaleron is therefore invariant under spatial translations and rotations. In this 
case, all zero modes must be of pure gauge origin, but the gauge is completely fixed. 

For n = 1, the vacuum operator Mq 00,1 has three zero modes. As is seen from Q5.14p , 
their eigenfunctions are just constants discussed in the previous section: = 5£t p /2. 
Since the spectra of M scal and M FP coincide (but these operators act in the different 
spaces), Mq P also has three constant zero modes, a p = t p /2, which have been discussed 
above. 

The next step is to pick up the physical modes. The partition functions Zj spec- 
ified by ( |4. 16 ) and ( |4.23| ) involve the ratios of the determinants, such that some of the 
eigenvalues cancel. Notice that O implies that Det'(M) = Det'(M scal )Det'(M vec ), 
whereas Det'(M scal ) = Det'(M FP ). The partition functions therefore are [BHj 



Z 1 = exp(-S E [A^})- 



a 



Det'(Mf al /a 2 ) 
2zJ\uj\^ Det>(Mr c /a 2 ) 



Z 



Ait 



Det'(M^ cal /a 2 ) 
a 3 p 3 \ Det'(M$ ec /sL 2 ) 



(5.18) 



Thus the physical oscillator modes constitute the part of the spectrum of M vec that 
remains after the subtraction of all eigenvalues of M scal . Since all eigenvalues of M scal 
are contained in the spectrum of M vec (besides the zero modes of Mq 6C ), such a subtrac- 
tion results in the change of the degeneracy factors of the eigenmodes of M vec . Using 
( |5.16| ), ( |5. 17[ ) , one obtains the eigenvalues and the degeneracies, {u 2 ,d}, of the physical 
oscillators: 



MT C /M{ cal 
also 



{-2,1}, {1,4}, {n 2 + (T 2 -3,2(n 2 -a 2 )}, n > 3, a = 0,1, 2; (5.19) 



M 



vec 




C /M s cat : {0,-3}, K,6(n 2 -1)}, n > 2. (5.20) 

Notice that o in ( |5.19| ) assumes three values which correspond to the three "colours" of 
the "gluon" , and the coefficient 2 in the degeneracy factor refers to its two polarizations, 
whereas in ( |5.20|) all (2 spin) x (3 isospin) gluon degrees of freedom are absorbed by the 
coefficient 6 in the degeneracy factor. Taking into account the temporal part of the 
fluctuation operators, — J^, whose eigenvalues are the Matsubara frequencies, 47r 2 / 2 //? 2 , 
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we obtain the following contribution of each physical oscillator, {u 2 , d}, into the partition 
function: 

AIM?-))""* 

For oj 2 < or uj 1 = the term with / = in this product should be omitted, since the 
corresponding negative or zero modes have already been taken into account. Notice that 
the negative degeneracy for the zero mode in ( |5.20|) arises simply because this mode is 
not in the spectrum of A^q £C and only in that of M. S Q al . 

We are now in a position to give the formal closed expressions for the partition 
functions. For the fluctuations around the sphaleron we obtain 



X 



X 

/=— oo cr=0,l,2 n=3 



n n n^(^+- 2 +^ 2 -3j , (5.22) 



and for the vacuum 



1 a 3 /? 3 - f 1 /47r 2 Z 2 \r 3 ~ ~ [ 1 /4vr 2 / 2 2 \l 3(n2_1) fK no , 

6 Evaluation of determinants and the transition rate 



We use the zeta function techniques (see [P0[|-|P7[, and references therein) to regularize 



and evaluate the infinite products entering Eqs.( [5.22j ), ( |5.23|) . The key steps of our 



analysis are presented in this section, whereas a large number of the technical details are 
given in the Appendix. 

The basic zeta function relation reads 

J] (jT = exp{-C'(0) - In// C(0)}, (6.1) 

where 

c(«) = Errb = fn / ^ 1 E^ex P (-a n )dt, (6.2) 

which in fact should be regarded as the definition of the product. 

We start by applying this to the zero mode contribution in Q5.23j ): 

The scale factor // is given by /i = sl(3/2tt, and the zeta function is 

00 1 

Ms) = -G y E~ = -6( R (s), Co(0) = 3, Co(0) = 31n27r, (6.4) 

1=1 1 
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where (r{s) is the Riemann zeta function. This gives 

-3 



n 



o2 



(6.5) 



such that the overall contribution of the zero modes together with their Matsubara 
excitations into ( |5.23| ) is 



3/33 oo 



a 3 /? 
Air 



n 3 



1 /47T 2 F 



/3 2 



1 

47T* 



(6.6) 



Next, consider the product in (5.22) due to the negative mode. The corresponding 
zeta function is 



C-(*) = £ * 2 



i=i 



2tt 2 



C-(0) = -|, C(0) = -ln 



'2y/2n p \ 
— sm(^)) (6-7) 



(see Appendix, Eq.( A.18f) ). The normalization factor is /i = (/3a/27r) 2 , which yields 

1 /4vr 2 / 2 V 



2 — n 



(3 2 



4 sin — = . 

V2 



(6.8) 



Now we want to take into account the contribution of the positive field modes. 
We introduce the spatial zeta function associated with the positive physical oscillators 

dm), O) 



(s P at(s) = 2+ £ 



n 2 — a 2 



<T=0,1,2 n=3 



(n 2 + a 2 - 3) f 



00 n 2 - 1 

-3E 



n=2 



such that 



One has 



eee 




Det'Mf al DetM 
\ DetM\ ec Det'Mf al 

1 



exp{CLt(°)} 



j Z" 00 i 

Cspat( s ) = — , r J (it t S ©spai(^) ? 



(6.9) 

(6.10) 
(6.11) 



where the heat kernel is 



e spat (t) = 2ex P {-t}+ E E^-^^pI-^+^-^I-sE^ 2 - 1 )^-^ 2 }- 

<t=0,1,2 n=3 n=2 

(6-12) 

We shall need the asymptotic expansion of this function for small t 



E C rt r - 



(6.13) 



( 4 ^) 3/2 ,=o,V2,i, 

The computation of the coefficients C r for r < 2 is performed in the Appendix: 

Co = C 1/2 = C 1 = 0, C 3/2 = -2(4vr) 3 / 2 , C 2 = 22vr 2 . (6.14) 
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Next, we introduce the thermal zeta function related to the spatial zeta function 

1 r°° 00 f /27r/\ 2 1 
Cfi(s) = f^J ^ E exp - (^— J t e spat (t)dt. (6.15) 

As a result, we can collect all parts together and represent the expression for the sphaleron 
transition rate in the following form: 



r = ^ = o /» 2 •//=»/ 6XP ^ + C ^ + llia C ^ ' 6 ' 16 

vr Z 8V27T 2 sm(/5/V2) [ <r ^ J 

Here the prefactor on the right hand side includes the contribution of the zero modes 
and the negative mode. The exponent contains the one-loop contribution of the positive 
oscillator modes, and the classical sphaleron action 

3tx 2 

S e [AM] = —f3. (6.17) 

g 

To compute the quantities Cg(0) an d C/?(0) entering Eq. ( |6.16|) , we represent the heat 
kernel in ( |6 . 1 2j ) symbolically as 



Os P at(t) = $>xp{-a; 2 t}, u 2 > 0. (6.18) 
Then the values of C/?(0) and Cg(0) are given by (see Appendix) 



167T 2 ' ^ V ' [ 167T 2 ^ V 2\ 

where the coefficient C 2 is specified by ( |6.14j ). 

All this allows us to represent the transition rate ( |6.16| ) as 



(6.19) 



' expl-^/?-^-^!-^))- (6.20) 



8V2tt 2 sin(/3/V2) \ 2 2 (a 
In this expression, the renormalized gauge coupling constant is 

?b = ?R"i^ ln ©- (6 ' 21) 

Here we have returned to the dimensionful a and replaced g by g(a ), where a = 1/yUo- 
This expression agrees with the renormalization group flow, such that it does not depend 
on the scale ao if g{sLo) is chosen to obey the Gell-Mann-Low equation. To fix the scale, 
we assume that the value of g(ao) is determined by the typical energy of the physical 
processes in the universe, that is by the physical temperature T(a ) = l//3a . Then we 
use the QCD data (see for example |p9| ) 

T(a ) = 100 GeV, = 0.12, (6.22) 

4"7T 

and assume that the weak coupling region extends up to some a max . One can choose 

^max ~ 10 100a . 
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Figure 1: The thermal function exp(— @(Fi — Fq)). 



So is the contribution of the zero field oscillations, that is, the Casimir energy, 

So = PP( spa t(~) + y (hi2 - 1). (6.23) 

It is worth noting that this quantity can be computed exactly in this case. The corre- 
sponding computation itself presents some methodological interest and is given in the 
Appendix. The result is 



£ = ^ + ^(ln2-l- 7 )+ f'dzVT 
6 4 Jo 



z 2 x 



x J tan (y) {{z 2 + A)F(z, t) + (Az 2 - 2)G(z, t, V2) + 9z 2 G(z, t,V3)}, (6.24) 
where 

„. . smh(Tczt) , sin(iiqzt) 2 sin(7rt) 

F(z, t)=t- . h ) G(z, t,q) = t- \ ; + - \ } 6.25 

smh{TTZ) sm{7iqz) 7T (1 — q z ) 

The numerical value is 

Sq = -1.084... (6.26) 
The contribution of the thermal degrees of freedom in ( |6.20| ) is 



/3(F 1 -F )=41n(l-e^) + 2 £ J> 2 - a 2 ) ln(l - e 



/Vn 2 +cr 2 -3^ 
o-=0,l,2 n=3 



-6^(n 2 -l)ln(l-e^ n ), (6.27) 

ra=2 

and the remaining sums in this expression can be evaluated numerically (see Fig.l). 

One can see that F% and F are precisely the free energies of the physical oscillators 
flOp , ( ^20|) . Altogether Eqs. (|6~20D - (|6T27|) provide the desired solut ion of the one-loop 
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log t0 (l//3) 



Figure 2: The sphaleron transition rate 



sphaleron transition problem. The numerical curves of T(l/f3) evaluated according to 
these formulas for several values of a are presented in Fig. 2. 

This solution makes sense under the following assumptions: 

1 1 3tt 2 . 
a < a„ < - < 6.28 

The first condition is the the weak coupling requirement. When the scale factor a is too 
large, the running coupling constant ( |3.21| ) becomes big (confinement phase), and the 



effects of the strong coupling can completely change the semiclassical picture. That is 
why our solution can be trusted only for small values of the size of the universe. The 
other condition in ( |6.28| ) requires that the thermal fluctuations are small compared to the 



classical sphaleron energy, such that the perturbation theory is valid. Notice that each 
curve in Fig. 2 develops a maximum at some temperature l//3 ma:r (a); however, this value 
seems to be already beyond the scope of the approximation: 3n 2 / ' g 1 : (a) ~ 2/f3 max (a). 
The subsequent decrease of T is presumably fictitious. Indeed, it is natural to expect 



that the transition rate is increasing with growing temperature |I2| . Thus, our results 
can be trusted at best only for 1/(3 < l//3 max (a). 

One can also find the high temperature limit for the solution (but the upper bound 
in ( |6.28j ) is to be assumed). To determine the asymptotic behavior of the free energy, the 
procedure is the following |3"Ifl . First, one returns to the zeta function (p(s) and replaces 
in ( |6.15| ) the heat kernel spa t(£) by its asymptotic expansion ( |6.13| ). Then one takes the 



integral over t, and the sum over I reduces to the Riemann zeta function. As a result, 
one arrives at the following asymptotic expansion for small /3 



^ (0) = 45^ + 2^^ + I^~i^ ln/5+ ^( 7 + ln 4^J 



+ 4^ E W ( £ ) U2r)T(r) + (^(0). (6.29) 
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This gives the free energy (see Eq.( p.l9| )) 



f = ~m + (1 "g 1 ;, 2 ^ 2 - ppc^c-^). (6.30) 



As an illustration we apply this to the vacuum fluctuations alone. The corresponding 
spatial zeta function is given by the last piece in Eq. (|6.9|) : Cs pa t( s ) = 3Cr(2s— 2)— 3(r(2s), 
its coefficients C r are computed in the Appendix (one has C2 = 0), and the result is 

Here 2tt 2 is the volume of S 3 (remember that j3 = 1/aT), and the leading term is just the 
free energy of a gas of noninteracting, massless particles with 2x3 polarization states. 

Now, let us return to the full expression ( |6.9| ) for ( S pat( s )- The values of C r in this 
case are given by fl6.14|) . Using ( |6.23|) , (|6.30|) one obtains 

- f3(F 1 - F ) = 41n/3 + C^(0) + ^(7 + ^)+ P + 0((3 2 ). (6.32) 
The quantity C spat (0) is computed in the Appendix: 

k = exp{Cl pat (0)} = ^ sinh 4 (vr)| sin( V2tt) | exp - I(y/2) - I(V3)} , (6.33) 

where 

rx rx / 2,t \ 

J(x) = ttJ t 2 coth(nt)dt, Z{x) = J f- — — + nt 2 cot(vrt)J dt, (6.34) 
with the numerical value 

k = 1250.21... (6.35) 

(It is interesting to observe the large value of k. The corresponding quantity for the 
electroweak sphaleron is suppressed by several orders of magnitude 0, ||.) 
Finally, one obtains (see Fig. 3) 

exp{-/3(F 1 -F )} = fi:/3 4 fl + ^/51n/5)+O(/5 5 ) as (3 -> 0. (6.36) 



As a result, we arrive at the following expression for the sphaleron transition rate 
in the high temperature limit: 

WH^expj-i^}, (6.37) 

where we have neglected also the Casimir term. As one can see from Fig.3, the high 
temperature approximation can be reasonable for 1//3 > 10 3 . On the other hand, the 
temperature should be less then the sphaleron energy 37r 2 /g 2 (a). These two conditions 
imply that ( |6.371 ) makes sense only for small (7(a): 10 3 < 37r 2 /g 2 (a), that is, for small a. 
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Figure 3: The ratio between the thermal function exp(— (3{F\ — Fq)) and its high tem- 
perature asymptote 1250. 21/3 4 . 

7 Concluding remarks 

In this paper we have obtained the exact solution of the sphaleron transition problem for a 
pure non-Abelian gauge field in a static Einstein universe. This has been achieved by the 
straightforward diagonalization of the one-loop fluctuation operators with the subsequent 
computation of the functional determinants in the zeta function regularization scheme. 
To carry out this program, the following points have been crucial: the high symmetry 
of the sphaleron solution and the fact that the sphaleron configuration consists of the 
gauge field alone. Actually, these properties of the model under consideration make it 
somewhat similar to the instanton theory M, [EH], J26||. It is worth noting that the 



solution obtained in this paper is unique in the sense that no other exact solutions of 
the sphaleron models in 3 + 1 spacetime dimensions are known. 

Eqs. (|6.20|) -( |6.27|) and ( |6.37|) are our principal results. They specify the number of 



transitions between the neighboring topological sectors per unit conformal time r\. It 
should be stressed that such transitions do not lead to any violation of chiral fermion 
number unless the thermal equilibrium between the different topological sectors is bro- 
ken. This can arise, for instance, when a fermion asymmetry is present. Then one has 
to introduce a small chemical potential fi for the fermions 0, ||. This favors those 
transitions which erase the asymmetry. Specifically, let AiV = Np — Np be the fermion 
number of the universe, then 

±AN^-t-T. (7.1) 

Notice that - since T defines the number of transitions in the whole universe - AiV refers 
to the whole universe as well, and not to the unit volume. For one doublet of chiral 
fermions, standard thermodynamics gives 

V r P* d P 3VT\ , « 2 , 2 , 
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where V = 27r 2 a 3 is the volume of 3-space, such that AiV = 27r 2 a/i/3/3 2 . Finally, passing 
to the physical time t = arj, one obtains the fermion number diffusion rate 

1 V = -^r(I). (7.3) 



AN dt 2vr 2 a y (3' 

In fact, we do not specify the nature of the fields under consideration. The dis- 
cussion of the possible applications of the results obtained in this paper will be given 
separately. At present we just mention where our results can be used. Let us recall the 
typical values of the parameters a, (3 and T = l/a/5. The range of T is restricted by the 
condition of the validity of the semiclassical picture: T > 1 4- 10 GeV ( "deconfinement 
phase"), whereas the metastability condition requires that the conformal temperature, 
1/(3, is not too high (see Eq.( |6.28[) ). The size of the universe, a = l/T/3, therefore should 
not be too large compared to 1/T. Such conditions can be met in the context of the 
finite volume QCD (the typical volume in that case is a 3 ~ 1 fm 3 ; see [[01 and references 



therein). Another natural possibility relates to the pre- inflation cosmology. In this case, 
the gravitating Yang-Mills field can arise in the context of a superstring theory |$8 1 . In 



fact, the semiclassical sphaleron transition picture applies also after inflation. However, 
the conformal temperature is enormously large then, 1/(3 = aT ~ s 1 / 3 ~ 10 28 , where s 
is the total entropy of the universe. Equivalently, one can say that the temperature T is 
huge in comparison with the sphaleron barrier 37r 2 /g 2 (a)a, such that there is no suppres- 
sion for transitions between the different topological sectors at all. Unfortunately, there 
are no reliable methods for computing the transition rate T in this limit. One can use 
a dimensional argumant to estimate that the rate related to the unit physical 4-volume, 
T/a 4 , should be proportional to T 4 . Then Y{1/ (3) ~ 1//3 4 , and (|7.3| ) yields 

' d AN ~ T. (7.4) 



AiV dt 



This agrees with the usual estimate for the fermion number dissipation rate at very high 
temperatures f§]. 
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Appendix. Zeta function techniques 

In Sections AD below the detailed computation of PP£ spai (— 1/2) and Cs pa t(0) ^ s P re ~ 
sented. The basic idea is to expand these quantities into certain series of the Riemann 
zeta functions and then to use the appropriate summation formulae. The asymptotic 
expansion for the heat kernel ® S pat{t) is derived in Sec.E. Section F contains the com- 
putation of C/3(0) an d C,a(0) f° r a thermal system. 
A. Summation formulae 



20 



Consider the generating function for the Bernoulli polynomials (see [57], p. 804) 

^T = E^(t)|, \x\<2k. (A.l) 



k=0 

Putting here t — one obtains 

1 00 (2x) 2fe 

coth(x) = - B ^^ttt, \x\ < 7T, (A.2) 
x k =o \ lK ) - 

where Bk(0) = B^ are the Bernoulli numbers. Using the relation between the Bernoulli 
numbers and the Riemann zeta function (see |)57| , p. 807), 



and considering the replacement x — > «x, one finds 

o oo / T \2k 00 / T \ 2k 

ooth(aO = -- EC- 1 )*!-) U2A:); cot(x) = — £ - Cb(2*)- (A.4) 
x fc f W Xk%^ J 

The sums on the right hand sides converge only for |x| < 7T, whereas the left hand sides 
are meramorphic functions on the whole complex plane. One can therefore consider these 
relations for any x as a result of analytic continuation. 

Let us now restrict ourselves to the real values of x. Integrating both sides of (]A.4|) 
one obtains 

00 (— ~\\k /r r \2k ™ oo I / v 2fe „ 

In the second of these formulae one should assume that x < n, unless the integration 
rule for the poles of cot(x) is specified (see below). 

Multiplying (|A.4j) by x 2 and integrating from zero to x one obtains 

|^g) 2 W2) = ^), 2) = ^(£), (A.6) 

where 

J(x)=tt t 2 coth(iTt)dt, (any x); J(x) = -n t 2 cot (irt)dt (x < 7r). (A. 7) 
Jo JO 

Eqs. (|A.5 )- (|A.7 ) will be used below. To proceed further, we return for a moment to 
Eq.( [A.l| ). Consider the equation which is obtained from (|A.1|) under the replacement 
x — > —x. Taking the difference of the two equations one obtains 



+ — = 2y4 + i(i)7T rr- (A. 

-1 e~ x -l + KJ (2k + l)\ 



k=0 



Using the explicit form of the k = term on the right hand side, 2B\(t) — 2t — 1, and 
utilizing also the following relation (see |57J, p. 807): 

( R (2k + 1) = (~l) fc+1 2 l J + ^ B»+i(*) cot(Trt)^, (A.9) 
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one finds 

°° / T \ 2k f i / xt -xt \ 

Replacing here x — > 2x and t — * 2£ — 1 , and considering also x — > ix one arrives at 

E(-i) fc (^) c*(2* + i) = -£/(*), E(r) U2fc + i) = ->), (A.11) 

fc=l 
where 

= ri / _ smh^\ ^ /ttA f = ,1 / _ sin^A ^ ^ ^ 
jo \ sinh(x) J V 2 / Jo \ sin(x) ' 



k=i 2 fc= iW 2 



Next, one deduces from (|A.11|) that 



7T 



£(_l)VV fc C«(2A; + 1) = --/(vr^). (A.13) 
fe=i 2 

Using 

3/2 



2 2 



(A.14) 



T(A; + 2) v^Fio T(fc + 3) 3v^Fio V 

we finally obtain the following formulae: 

g(-l)V fc ^±^ C fl (2fe + l) = -2v^F f dzy/T=*I(*uz), 
kTi T ( k + 2) Jo 

E(-l) fc ^ r r( fc + + 1 3) ) C fi (2fc + 1) = -1 V5F jf 1 (l - , 2 ) 3/2 /(tti/z), (A.15) 

together with the corresponding two relations obtained by v — > iv, I(ttvz) — > I(irvz). 
Here the integrals can be considered for arbitrary values of u, which defines the analytic 
extension of the series. 

B. The basic zeta function relation 

Consider the following zeta function 

C°(*) = E , 2M tTs = / t 3 " 1 E ex P {-^ 2 + v 2 )}dt. (A.16) 

^ (n 2 + z/ 2 )* r(s) Jo J 

In order to express this function in terms of the Riemann zeta function, we expand 
exp(— tv 2 ) = J2k{~~ ]-) k v 2k t k /k\, and perform the integration over t. Then the sum over 
n gives the Riemann zeta function, such that 

C(s) = ( R (2s) + £ ^^ fc£ T(^^(2fc + 2s). (A.17) 

This relation will turn out to be especially handy. Taking the pole of the gamma function 
at s = into account, using Eq. flA.5p and utilizing (r(0) = — |, Cr(0) = — |hi27r, one 
finds 

tf(0) = -i, s«°' = - h ^ + S^«^ = h 2S£(S0- (A ' 18) 
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For a harmonical oscillator, for instance, one has 

(A.19) 

which gives rise to the formula ( |6.8|) in the main text when uj = iy/2. 
We shall also use another zeta function 

CM - t j-^y, - W - 2) + 1 tf^J±+A U2k + 2s _ 2) . (A . 20) 

Using Eqs. (|A.6| ) and (|A.7|) one obtains 

C(0) = 0, ^C 2 (0) = 2&(-2) + ^(i/), (A.21) 

where J7"(z/) is defined by ( |A.7| ). 

C. Evaluation of PPC P at(-|) 

We consider the spatial zeta function 

00 n 2 - a 2 00 n 2 - 1 

C P at(s) = 2+ £ £ - -3^^-, (A.22) 

which can be represented in the form 

Cs P at(s)=A(s) + B(s), (A.23) 

where 



A(s) = 2 + 3Ca(2s) - 3( R (2s - 2) + J> 2 - 1) (4 + v 2 ) , (A.24) 

and 

= E {C°(* - 1) - (1 + - E(3 + 2^ 2 ) {£(«) - (1 + v 2 )- 3 } ■ (A.25) 

Here C°( s ) is defined by Eq. (|A.16|) , and v 2 = a 2 — 3 = 1, —2, —3. In A(s) one can simply 
put s = -1/2, which yields PP.4(-l/2) = ^4(— 1/2) = -91/40 - 3y/2. 

Let us now consider B(s — 1/2) in the limit s — > 0. We first analyze — 1/2). 
Using Eq.( |A.17|) we obtain 

= C R (2,-1)-, 2 ^|yCk(2 S +1) + £ ^u 2 ^ k + S -^ 2) U2k + 2s-l). 

(A.26) 

Let us now consider the limit where s tends to zero. Then the second term on the right 
hand side diverges due to the pole of the Riemann zeta function, 

(r(1 + s) = ± + 1 + 0(s 2 ), (A.27) 

where 7 is the Euler constant. The remaining terms in ( |A.26| ) are all finite. The principal 
part is 

pp ^ML"u(<°0)- < A - 28) 
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Using Cr(~ 1) = — 1 /12, r(— 1/2) = and replacing k — > k + 1 in the sum entering 

Eq.( |A.26|) , one obtains 

PPC°(--) = -— - -u 2 (l - 7 ) + -^-y(-l) k u 2k T{k ^ 1/2) ( R (2k + 1). (A.29) 
SlA 2 J 12 2 v 2A^i r(fc + 2) sm ; v ; 

Finally, taking into account Eqs.( A7T5| ), one arrives at 

PPC°(-i) = ~ - ^ 2 (1 - 7) - ^ 2 jT 1 ckvT^JCTri/z), (A.30) 
where I(x) is defined by Eq.( [A.12| ). Similarly, one obtains 

u 0/ ^\ 1 1,2 1 ,4/i 3 \ 4 Z' 1 , /-, 2\ 3 / 2 



PPC U (--) = — - -z^ - -zAl - -7) - z/ 4 / dz(l-z 2 )' IOkvz). (A.31) 
SM 2 ; 120 8 2 v 4 ' ; Jo v J K J K J 

The next step is to insert these relations in Eq. (|A.25 ) and to compute the sum over 
v 2 = 1, —2, —3. The v 2 = 1 case presents no problems, whereas the negative values of v 2 
should be treated with some care. Let us pass in ( |A.30|) to negative z/ 2 , v — > iq, where 
q = \q\. For q < 1 one obtains 

PPC°(-i) - PPC° (~) = -^ + \q 2 (l-l) + q 2 [ dz>/T=#I(icqz), (A.32) 

where I(x) is defined by Eq. ([A.12 ). Next we consider the following integral representation 
for y/1 — q' 2 which is valid for < q < 1: 



7 =1 _??!/ 1 (fev /T^/ 1 i^ tm (^W (A.33) 
n Jo Jo 1 — q 2 z 2 V 2 / 



This implies 



1 1 



PPC°(-l/2) -v /l -^ = - 1 -TVJ + 2 ?2(1 ~ 7) + 



+ q* /* = * f ( t - ^4 + 2 -^^) tan (?*) dt. (A.34) 
jo Jo \ sm{7Tzq) tx 1 — q A z A J \ A / 

Here we can safely extend the range of g from 0<g<lto0<g<2, thus taking the 
values q = \/2,\/3 into account, that is, v 2 = —2,-3. This allows us to compute the 
second sum over v 2 in the formula ( |A.25| ) for PPB(— 1/2). The first sum can be done 
with a similar rearrangement, using the formula 



Finally, collecting everything together, we arrive at 

PP(s P at(-h = ^-^7 + 5^ 1 (l)-^ 3 (l) + +2P 1 (v / 2)-4P 3 (v / 2)+9P 1 (v / 3)-9P3(v / 3); 
2 o 4 

(A.36) 



where 



K„M =/>(!- -- 2 )"" 2 / (« - Sgf) tan (£) *, 



^ / n Z" 1 , 9 \m/2 /■! / sin(7T2;(7t) 2 sin(7rt) \ /7rt\ , , . 

= j[ * (i - , 2 ) / (* - -4^ + - I r3 7 4) tan (t) *• (A - 37) 
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D. Evaluation of (' spat (0) 

The procedure in this case is similar to that of the previous section. One again starts 
from Eqs.( [A.2"3"|) -( [A.25"D , but now, using (|A.16|) and ( |A.20| ), it is convenient to represent 
the function B(s) in the following equivalent form: 

Bit) = e {C(s) - (i + v 2 r s ) - E(s + » 2 ) {C(s) - (i + v 2 r s ) . (a.ss) 

ifi 

The next steps are straightforward. Using ( |A.18|) and (|A.21|) one obtains 

v{\ + v 2 ) 



B'(0) = E {2CA-2) + J{y) + ln(l + v 2 )} - E(3 + v 2 ) In 



2 sinh(7rz/) ' 



(A.39) 



where J7"(V) is defined by ( |A.7| ), and v 2 = 1, —2, —3. For negative u 2 = —q 2 < 



one 



again needs some minor modifications. No changes are needed for the second sum in 

q{l-q 2 ) 



In 



In 



2sinh(7rz/) 2sin(7rg) 
where we can put q = y/2, In the first sum one has 



J(u)+\n(l+u 2 ) 
where 



J(q)+\n(l-q 2 ) = -n t 2 cot (nt)dt- 



i 2tdt 
o 1-t 



2t 



+ Tit 2 COt(7Tt) J dt 



lo VI - 1 2 

which can also be extended to the values q = ^/2, v 7 ^- For *A'(0) one has 



(A.40) 

2 _ -X{q), (A.41) 
(A.42) 



^'(0) = -3 ln(27r) - 6^-2) + 3 In 2, 



(A.43) 



which finally gives 



CU(°) = ln (^T sinh V)l sin(v^r)|) + ,7(1) - J(v^) - T(yfi), (AAA) 

with J7"(a;) and T(x) being defined in Eqs. (|A.7|) and ( A.42 ), respectively. 

E. Asymptotic expansion of Q spa t(t) 
Let us consider the following function 



e(z,j,^) = £(n 2 -j)e- 



-(n 2 +k)t _ _ e ~kt 



d \ » . 



(A.45) 



We want to find its asymptotic expansion for small t 

1 



e(i,j,k\t) 



(4vrt) 3 / 2 



E w 

r=0, 1/2,1,... 



(A.46) 



Such an expansion can be obtained with the use of the theta function identity j36] 



oo oo f ^2 

E exp(-tn 2 ) = W - E expj-— n 2 

n=— oo n=— oo I 



(A.47) 
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which implies that 

]Texp(-tn 2 ) ~ -J exp(-tn 2 ), (A.48) 

n=i n=l 

since for t — > all terms with exp{— 7i 2 n 2 /t}, n ^ vanish faster than any power of i 
and can therefore be omitted. Inserting this into ( |A.45|) and expanding the remaining 
exponents we find 

©(<, J, ~ ^ - + 2j) +j/2 + E (J - n 2 ) + ^k(k + + 0(t). (A.49) 



The coefficients C r for r < 2 are therefore 

n=l 

(A.50) 



Co = 2vr 2 , C 1/2 = 0, Ci = -2vr 2 (A;+2j), C 3/2 = Stt^^ + E (j-n 2 )), C 2 = 7T 2 k(k+Aj). 

1 n=l 



The heat kernel (|6.12 ) in the main text can be represented as 

e spat (t) = 2e-* + E e(3 )( x 2 ,(T 2 -3|t)-3e(2,i,0|t), (A.51) 

(7=0,1,2 

which finally gives rise to Eq. (|6.14)) . 
F. Evaluation of £s(0) and ££(0) ||35l 

Consider the thermal zeta function for an arbitrary system of harmonic oscillators with 
positive energies 

1 r°° 00 f /2-7r/\ 2 1 

U*) = f£) / ^ z £ ex P - (— j * | (A-52) 

where 

6(t) = £exp{-u; 2 t}, w 2 >0. (A.53) 

Transforming the sum over / in ( |A.52| ) with the use of the theta function identity (|A.47 ) 
and using the integral representation for the Kelvin functions 

^w-KDT*^ 1 ""!-*-^} 1 (A - 54) 



one obtains |}5] 

= 2^FM y -< S - 5> + A I? (Q"' 2 ^/-^), (A.55) 

where 

npat(s) = / dt f^Qit) = Cspat(s)T(s). (A.56) 
This function has the following pole structure |3c 



« = (SF?JT^2 +/(!) ' (A ' 57) 
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where C r are defined by the asymptotic expansion of O(t), and f(r) is an entire analytic 
function of s. This relation implies 

Y spat (s - \) = + PPY^i- 1 -). (A.58) 



Taking Eq. (|A.56 ) and the properties of the gamma function into account, 

-L = s + 7S 2 + 0(s 3 ); T{~ + s) = -2v/?{l + (-7 + 2 - 21n2) S } + 0(s 3 ), (A.59) 
one therefore obtains 

1 Y spat (s -h = £i + { (1 ~ c ln ?° 2 - PPC**(4)1 * + 0{s\ (A.60) 



2^T(s) bvmy 2 J 16vr 2 ' { 8vr 2 2 
Finally, using 

nr 00 i 

K 1/2 (z) = ^-e-\ J2 -e~ l * = - ln(l - e~% (A.61) 

one arrives at 



16tt 2 (V 8tt 2 2 / 

which gives rise to the formula ( |6.19|) in the main text 



(A.62) 
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